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We investigate the nonlinear Meissner effect microscopically. Previous studies did not consider 
a certain type of interaction effect on the nonlinear phenomena. The scattering amplitude 
barely appears without being renormalized into the Fermi-liquid parameter. With this effect 
we can solve the outstanding issues (the quantitative problem, the temperature and angle 
dependences). The quantitative calculation is performed with use of the fluctuation-exchange 
approximation on the Hubbard model. It is also shown that the perturbation expansion on the 
supercurrent by the vector potential converges owing to the nonlocal effect. 

KEYWORDS: nonlinear Meissner effect, unconventional superconductor, electron correlation, spin fluctu- 
ation, vertex correction, nonlocal effect 



Many unconventional superconductors have recently 
been discovered. The evidence that these superconduc- 
tors are non-s-wave is obtained by thermodynamic mea- 
surements in most cases. These types of measurement 
provide only information averaged over the Fermi sur- 
face and then other measurements are needed to deter- 
mine the position of nodes. Detailed information about 
the node is useful for judging the accuracy of various 
theories. One of these, the nonlinear Meissner effect 
(NLME) , which provides a measurement of the magnetic 
field (if) dependence of the magnetic field penetration 
depth (A + SX(H)), was proposed by Yip and coworkers 
(YS). 1 ' 2 Their proposal is based on the Doppler-shifted 
energy spectrum and its predictions are summarized as 
follows, (i) The supercurrent has a nonanalytical form 
(written as A\A\, A is the vector potential) and SX(H) 
is proportional to the magnetic field \H\ as a result, (ii) 
5X{H) varies with the direction of the applied magnetic 
field and therefore SX(H) provides information on the 
position of nodes. 

To date, the experiments have not provided deci- 
sive results because the effect is small and tends to be 
masked by many extrinsic effects. The first investigation 
of NLME was carried out by Maeda et al., 3 however this 
experiment was performed with the magnetic field per- 
pendicular to the CuO plane. The precision was also poor 
(order of 10A) and, at present, the observed quantity is 
considered to reflect extrinsic effects. Experiments with 
high precision (order of 0.1 A) were carried out by Bidi- 
nosti et a/. 4 and Carrington et al. 5 The most reliable 
results in ref.4 are summarized as SX(H) oc if 2 , the tem- 
perature dependence is weak and the angle dependence 
is not observed. These contradict YS's theory. YS also 
predict that 5X(H) oc if 2 below the crossover field, how- 
ever, the temperature dependence is strong. Two groups 
have attempted to detect the transverse magnetization, 
which is the other prediction made by YS. 6 ' 7 The recent 
experiment with higher precision (two orders of magni- 
tude) 6 showed that the amplitude of this quantity is at 
most one third of the predicted one and it is almost at the 
measurable limit. Extensions of YS's theory were made 
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by several researchers. 8-11 The common results obtained 
by them are that the theoretical predictions are incon- 
sistent with the experimental results. Therefore, some 
papers suggest that the experiments observed extrinsic 
effects. (In ref.5 the possible extrinsic effects are listed 
(the vortex contribution, the weak links, the interlayer 
contribution).) The important point is that the values 
predicted by YS and other researchers are at least of 
the same order of magnitude, compared with the values 
obtained in the experiments. Therefore, it is not possi- 
ble that the intrinsic theoretical value is masked by the 
extrinsic effects which are suppressed up to 0.2A. 

Then, the question arises as to whether the exist- 
ing theories are correct for judging experimental results. 
Here, we discuss the NLME effect on the basis of the 
perturbation theory and show that the previous theories 
have some defects. The perturbation expansion by the 
vector potential on the supercurrent and the magnetic 
field penetration depth converges owing to the nonlocal 
effect. The intermediate-states interaction 12 (electron- 
electron) which is not included in the conventional qua- 
siclassical approximation 13 exists and makes a domi- 
nant contribution. This effect solves the inconsistency 
between the theory and the experiments on the value 
of SX and its angle and temperature dependences. We 
adopt the fluctuation-exchange (FLEX) approximation 
for the quantitative calculation. The many-body effect 
on the response function is included on the basis of the 
conserving approximation. 

Our theory is based on the evaluation of the response 
function in the supercurrent which is exactly expanded 
by A up to the third order. The expression for the su- 
percurrent is written as 

J,{q) = -K$(q)A v (q) - f K$ a (q,q')A a (q')A v (q - q') 



41/3(9- <?', q")A (q")A a (q> - q")A v {q - q'), (1) 



q',q" 



where iff 1 ' 2 * 3 ) are the response functions in the pertur- 
bation expansion, fj,, v, ... are the spatial dimensions and 
the summation of the repeated indices is taken. ap- 
pears in the usual linear response theory. The K^ 1 term 
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is dominant in the magnetic field dependence of A be- 
cause vanishes. 

First, we show the convergence of the perturbation ex- 
pansion. By analyzing various terms it is shown that 
the most divergent term in the local limit is the type 
(a) term in Fig. 1 and is written as, K^ 3a \q,q',q") = 

V Sfe.n V k Tl l G k+q/2Gk+qi-q/2Gk+q"-q/2Gk-q/2}, where 

T, Vk and Gk are the temperature, the velocity and 
Green's function in the superconducting state, respec- 
tively, and k = (fc,ie„) (e„ = 7rT(2n + 1), n is integer) 
and q = (q, 0). The uniform component is written as 
(0,0,0) oc -f-i (A is the maximum of the su- 
perconducting gap). This is the result of YS's theory. 
K^\q,q',q") diverges for T -» at q = q' = q" = 0, 
however, this term is integrated by q' and q" in the 
expression of the supercurrent and by q, q' and q" in 
the case of SX. Then, the behavior of K^'(q,q' ,q") 
in q-space comes into question. The q-dependence of 
K^ a \q,0,0) for small q and at T = is written as 
K^ 3a \q, 0, 0) oc — -j^ v \, (v± is the mean value of the 
the interlayer velocity with magnetic field parallel to 
the a6-planc.) It is difficult to determine analytically 
the dependence of K^ 3 \q,q' ,q") on q, q' and q" , how- 
ever the nonzero values of q' and q" do not make this 
term more divergent than 1/q and if we consider the 
symmetry relation, K^ 3 \q, q, q) = K^ 3 \q, 0, 0), the form 
for small q, q' , q" is considered to be, K^ 3 >(q,q' ,q") oc 

l/^d ) 2 + («' - f ) 2 + («" - §) 2 + (f ) 2 - The q- 1 di- 
vergence of K^ 3 \q,q' ,q") guarantees the convergence of 
the perturbation expansion on the physical quantities. 
Therefore, the nonanalytical behavior does not appear 
contrary to the prediction of YS. 14 

Next, we estimate the various terms in K<- 3 K (The di- 
agrammatic representation of is given in Fig. 1.) 
There are many cumbersome terms unlike the linear re- 
sponse case. 15 The approximation used is as follows. The 
three-point vertex correction connected with the odd or- 
der of the external field is omitted. This is because the 
velocity is an odd function in wave-number space and 
the integral is small. The same approximation holds in 
some types of the six- and eight-point irreducible ver- 
tices. If we consider a system with a strong momentum 
dependence such as the underdoped region and its dop- 
ing dependence, this type of vertex is necessary. 16 We 
consider, however, mainly the temperature and the an- 
gle dependence, and then the above terms have a slight 
influence. 

The formalism for deriving 5X(H) consists of the 
Maxwell equation with the specular boundary condition 
(here, we consider the situation where the external mag- 
netic field is applied parallel to the y-axis and the super- 
conductor occupies z > 0) d = 2HS(z) — ^-J x (z), 
(H is the applied external field) and the nonlinear Pip- 
pard equation (eq. (1)) From these two equations the 
nonlinear equation for A is obtained and is solved by the 
perturbation method (not self-consistently) because we 
use the perturbation method to obtain the response ker- 
nel K( 3 \ Then, the nonlinear correction for A is obtained 
as (the magnetic field is parallel to the intralayer crystal 




Fig. 1. Diagrammatic representation of K^ :i K The solid lines ex- 
press the propagator of the electron in the Nambu representation. 
The wavy lines express the electro magnetic fields. The vertex 
with n-wavy lines denotes the n-times derivative of the disper- 
sion of electrons and the circle at the vertex means a diagonal 
matrix T3 with Tr-?3 = 0. The vertex with a shaded triangle sat- 
isfies the integral equation with a four-point irreducible vertex 
/( 2 ). The shaded rectangle denotes the reducible four-point ver- 
tex r^ 2 ' . 1^ and /( 4 ) represent the irreducible six- and eight- 
point vertices, respectively. 



axis) 

xD xx {q")D xx (q' - q")D xx {q - q')D xx {q)(2) 

where the definition of the magnetic field penetration 
depth is A := i J °° H y (z)dz = -± J %LA x (q), and 

D xx (q) := l/(q 2 + ^K xx \q)). 

First, we consider the terms without the intermediate- 
states interaction. We classify these terms into two 
groups. One group consists of divergent terms without 
nonlocality and the other of terms similar to (the 
nonlocality is negligible). K^ 3a ^ and K^ 3h,c,d ^ are cat- 
egorized into the former and the latter groups, respec- 
tively. We estimate these terms only with the experi- 
mentally observed values. (On the other hand, a spe- 
cific microscopic model and approximation arc needed 
in the case of the intermediate-states interaction.) It is 
difficult to analytically calculate the paramagnetic term 
if( 3a ) at a finite temperature, however, by noting that 
Ek+q — Ek — vj_q near nodes and the characteristic value 
of q is A -1 it can be shown that there is a crossover 
temperature To ~ ~ (£ c is the interlayer 

coherence length). By substituting K^ 3a \q,q',q") into 

cq. (2), S\ (3a) = (-g;) \K (3a H^). Here, « (3a) (^) is 
a dimensionless quantity and its integration is calcu- 
lated numerically. (H c := cf) /2^/2T:^ ab \ is the thermo- 
dynamic critical field, </>o = 2irhc/2e, £ a b is the in-plane 
coherence length). K^ 3a ^(^) is strongly dependent on 
temperature and shows a maximum around T ~ To. If 
we put A = 1600A, T /Ao = 0.01 and H c = 8000G, 
<5A( a ' < 0.2A for H = 200G; this is small compared with 
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the experimental result <5A| cxp 
than one tenth of SX due to the intermediate-states in- 
teraction discussed below and is negligible. 

Next, we consider the intermediate-states interaction. 
(The local approximation holds as in the case of (b,c,d,e) 
and K^\) The terms are calculated after they are trans- 
formed into equations which explicitly equal when 
A = 0. The term with ((g)-type) is evaluated as 
follows. The reducible four-point vertex satisfies 

f ( 2 ) (jfe, k') = J( 2 ) (k, k') + ^J2 ii2) ( fc ' k "yg(k")f ( 2 ) (fc", k') 

k" 

(3) 

with the irreducible four-point vertex I^ 2 \k,k'). ( de- 
notes the 4x4 matrix in the particle-hole space.) We 
use eq. (24) in ref.17 as g(k). There are various terms in 
l( 2 \k,k'), the particle-hole vertex, the particle-particle 
vertex, the number-nonconserving vertices. These are 
the functional derivatives of the self-energy by Green's 
function (the conserving approximation 18 ). One of these 
vertices (the number-conserving particle-hole vertex) is 
written as follows with the FLEX approximation (the 
FLEX approximation in the superconducting state, for 
example, see refs.19 and 20) 

I c (k, k') = V£_ k , - Gk- q W q (G k >- q + G« +q ). (4) 
i 

Here, 



5A. 6\( h,c,d ' e ^ is smaller experimental value.) On the other hand, the experimen- 



yn 

g 



U 



3 {U X s q ? 1 (U X c q f 



21-Ux s a 2 1 + Ux c a 



(5) 



with x s q ' c = -yY, k ( G k+q G k ± F k+q F k ), ((s,c) corre- 
spond to (+, — ), and 



W„ = U 2 



+ 



2 (1 - Ux s q ) 2 2 (1 + U X c q ) 2 



1 



(6) 



The other vertices are derived in the same way. We adopt 
the Hubbard model with the on-site Coulomb interaction 
U and take the same dispersion of electrons as in ref.16. 
The terms with 1^ and 1^ are calculated in the same 
way (e.g., I^{k u k 2 , k 3 ) = 5 2 t k J SG k ,8G k2 ). We cal- 
culate without integrating out the incoherent part 
to derive a low-energy expression (e.g., ref.16.) There- 
fore, the effect of the renormalization factor is implicitly 
included. 

To compare the calculation with the experimental re- 
sults quantitatively, we consider the following quantity 



1 5\{H) 



iM 



a 2 K^ 



H 2 



X 



16 {e/hcfK^t/a ifM 



(7) 



Here, a is the lattice constant. We numerically calculate 
and by putting t = 1 and a = 1 and quantify 
the values of t and a in units of [eV] and [A] , respectively, 
J-^MM1 = _ .i3 x 10- 13 - 



\G- 



If 



1 SMH) 
X 



and then 

we put t — 0.25[eV] and a = 8 [A], we get 
0.35 x 10" 7 [G- 2 ] for U = 6.0 and the hole doping 
S = 0.20. (K^ ~ -12.8, ~ 0.1 and the dominant 
contribution comes from K {z 9). This value of yields 
A ~ 2600A. This is roughly 1.6 times longer than the 



tal result is 



1 S\(H) 



(0.7 - 1.0) x 10" 7 [G- 



21 4 



Our 



H 2 A 

calculation is quantitatively consistent with the experi- 
mental results in order of magnitude. As for the param- 
eter dependence, the value of jp is not strongly 
dependent on the parameters U and S in the FLEX cal- 
culation. For example, -^ S -^p- ^ 0.34 x 10- 7 [G~ 2 ] for 
U = 7.0 and S = 0.20 and 0.33 x lO^G" 2 ] for U = 6.0 
and S = 0.15. This is because the effect of the renormal- 
ization factor on and (K^) 2 cancels each other (as 
for [/-dependence) and the integral-equation structure 
for r^ 2 ) weakens the variation of the spin-fluctuation ef- 
fect on /( 2 ). If we put W q — > U 2 , this corresponds to 
the case of the weak spin fluctuation, for example, the 
more overdoped region, this results in a smaller value of 
_i^£A(^2 therefore, an experimental study on the dop- 
ing dependence is expected. 

To investigate the angle dependence we consider the 
case where the applied field is parallel to the node 
direction (6X450 (if)). In this CciSG, Kxxxx in cq.(2) is 
replaced by (K X 3 J XX + 3K x 3 J yy )/2. Then the relation- 
ship between K^^^ and K^ aa with ji 7^ a plays 
an important role in the angle dependence. If we con- 
sider a conventional s-wave superconductor, the rela- 



tion K 



(3) 



3K, 



(3) 



holds because < 



>FS = 



3 < vfot >fs in the superconductor with no nodes and 
= and v 2 ,, = 



v. 



nnn — v ™ — v^f.Vaa hold in the electron gas. 
(< ... >fs denotes the average over the Fermi surface 
and u w = dv^/dkfj,, etc.) Both of these relations do 
not hold in the unconventional superconductor in the 
lattice system and therefore the relationship between 
and K$ aa \^ a is not trivial. In fact, K X 3 XXX and 
K x 3 x a yy give the same contribution to 5X because a dom- 
inant contribution to the integral over the Fermi sur- 
face comes from nodes (v x = v y at this point) except 
for T » T . Then, <5A 45 °(ff) = 25X ab (H) in the con- 
ventional quasi-classical approximation. 21 We made sure 
above, however, that the intermediate-states interaction 
term contributes to SX sufficiently and can be dominant. 



In this case, K x 3 J yy = K xxxx /3 with the approxima- 
tion noted above and then SX 45 o (if) = 5X a b(H). This 
explains the experimental results. (The reason for this 
is that the correlation between different vertices is bro- 
ken by the intermediate-states interaction. For example, 

SqiSk GkV kil G k V kv G k G k - q W q / fe , Gfc' -qG k > V k > a G k > V k >/3G k > 

is negligible for /1 7^ v or a 7^ (3. The same holds for the 
case of I ( 3 ) and I^K This discussion also applies to the 
temperature dependence.) 

The temperature dependence of <5A is as follows. SX cx 
1/T for T > T in the conventional quasi-classical ap- 
proximation. On the other hand, the temperature de- 
pendence of the intermediate-states interaction term is 
same as that of (T-linear) and the decreasing rate 
compared with the value at T — is almost same. There- 
fore, 6X shows a slight increase. (If A increases 5A per IK 
as the experimental results indicate, the increasing rate 
of SX is ~ 0.015 [A/K]. This value is no larger than 
the experimental precision.) 

The other phenomenon related to the NLME is the 



(3) 
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transverse magnetization. The predicted behavior in 
rcfs.l and 2 is that the supercurrent is not perpendic- 
ular to H except for the case in which H is parallel to 
the nodes or the antinodes, and therefore the transverse 
magnetization has a period of it/ 2 as the direction of H 
is rotated. Our perturbation theory shows that the su- 
percurrent is written as follows in the arbitrary direction 
ofH. 

J,(q)^{K^(q)A(q)+ [ [ [K^(q,q',q")X 2 

J q' J q" 

+^% y (q, q', q")Y e 2 }A(q")A(q> q")A(q q')}X e . 

(8) 

Here, (Xg 7 Yg) = (cos#, sin#), (sin#, cos#) for ^ = x,y, 
respectively. 6 is the angle between the applied field 
and the intralayer crystal axis.) If the relation = 
holds, the transverse magnetization does 
not appear. Then, we can have the same discussion as 
in the case of SX. 

Finally, we comment on previous studies. The nonlocal 
effect considered in ref.9 is different from our approach 
in several points. The behavior <5A oc H 2 at low H is 
seemingly the same as that in the perturbation approach. 
They predict, however, SX oc H above the crossover field 
H * and argue that the NLME is unobservable owing to 
H* > H c i . They consider that the nonanalytical current 
exists above H * . Although they do not consider the angle 
dependence of SX and the transverse magnetization, their 
theory contradicts the experimental results. The origin 
of their error is that they consider K^{q, A q=0 ). They 
derive H* by comparing the effect of q and A, however, 
it does not make sense to compare the intrinsic spatial 
variation with the external field. The absence (or very 
small value) of the transverse magnetization below the 
first vortex penetration 22 implies the absence of H* . 

The quasi-classical approach in ref.23 gives observ- 
able values {SX ~ lA for H ~ 200 [G]) with the 
experimental parameters (H* ~ 2[T] and A/£ ~ 
100). Then, this also contradicts the experimental 
results in the angle and temperature dependences 
qualitatively. The cause is as follows. The interac- 
tion with the external field in the Gor'kov equa- 
tion ^ (V r + ^ - i^A(R + r/2)) 2 G(r, R) is approxi- 
mated as (V r + - i^A(R)) G(r, R) in the quasi- 
classical approach. (The propagator is transformed as 
G(x, x 1 ) -» G(x - x', z±f) = G(r, R).) This means that 
the external field interacts with the center of mass of 
the electron propagator and therefore the nonlocal effect 
is underestimated. The comparison with our (a)-term is 
as follows. The Green function in the third order of the 
external field is written as 

G {3) (k,q) -> G k+q/2 vA q2 G k+q/2 -q 2 vA qi - q2 

91,92 

~xG k+q /2- qi vA q - qi G k - q/2 (9) 
for the case of the nonlocal effect included correctly, but 
G (3) (fc,g) -» G k+q/2vA q2 G k+q/ 2- q2 / 2 vA qi - q2 

91,92 



xG k+q /2- qi /2vA q - qi G k (10) 

in the quasi-classical approximation. This is interpreted 
as meaning that the magnetic field penetration depth 
effectively doubles in this approximation and then SX 
roughly increases eightfold. Therefore, the quasi-classical 
term ((a)-tcrm in our paper) makes less contribution to 
SX if it is evaluated properly. 

Larkin and Ovchinnikov suggest that the quasi- 
classical approximation does not give correct results in 
some cases. 24 Our theory presents a definite example of 
this proposition. 

In this paper, we present the microscopic formulation 
of the nonlinear Meissner effect. We show that the pre- 
vious studies on this effect are insufficient and some of 
them are incorrect. The nonanalytical response is intrin- 
sically absent. The experimental results possibly observe 
the intrinsic NLME. This is not YS's one, but originates 
from the intermediate-states interaction. We consider 
that this effect is interesting because it does not appear 
in the zeroth order of interactions but it reflects interac- 
tions between quasiparticles themselves. The spin fluctu- 
ation is quantitatively dominant in our calculation. This 
is consistent with the properties of the high-T c cuprates. 
Experiments on other materials and the theoretical inves- 
tigations of various scattering mechanisms are expected 
in the future. 

Numerical computation in this work was carried out 
at the Yukawa Institute Computer Facility. 
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